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ABSTRACT 

We present the analysis of an extensive set of high-resolution spectroscopic obser¬ 
vations of HD 93161, a visual binary with a separation of 2". HD93161A is a spec¬ 
troscopic binary, with both components clearly detected throughout the orbit. The 
primary star is most probably of spectral type 08V, while the secondary is likely an 
09V. We obtain the first orbital solution for this system, characterized by a period 
of 8.566±0.004 days. The minimum masses of the primary and secondary stars are 
22.2±0.6 Mq and 17.0±0.4 M 0 respectively. These values are quite large, suggesting 
a high inclination of the orbit. The second object, HD93161B, displays an 06.5V((f)) 
spectral type and is thus slightly hotter than its neighbour. This star is at first sight 
single but presents radial velocity variations. 

We finally study HD 93161 in the X-ray domain. No significant variability is de¬ 
tected. The X-ray spectrum is well described by a 2T model with kTj ~ 0.3 keV 
and kT 2 ~ 0.7 keV. The X-ray luminosity is rather moderate, without any large emis¬ 
sion excess imputable to a wind interaction. 

Key words: stars: individual: HD 93161 (A,B) - stars: binaries: spectroscopic - 
stars: early-type - stars: fundamental parameters - X-rays: stars - X-rays: individual: 
HD 93161 


1 INTRODUCTION 

A few years ago began an extensive campaign to study the 
Carina Nebula, one of the largest concentrations of young 
and massive stars known in our G alaxy. This campaign 
was initiated by the X-Mega group llCorcoran et all 11999 0 
in order to gain a good knowledge of the cluster’s binary 
population prior to its observation in the X-ray domain by 
the ROSAT satellite. Such multiwavelength observations 
of spectroscopic binaries in open stellar clusters enable to 
constrain the fundamental parameters of the early-type 
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servatories (NOAO). NO AO is operated by the Association of 
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stars and to investigate the interactions between their 
components in more depth. 

The efforts were focused more specifically on the 


massive stars of the Trun 

ipler 16 (Trl6) open clus- 

ter (Albacete Colombo et al. 

20011: [Morrell et alJ 12001!: 

Pa.nw et al. 

l200ll: 1 Albacete C( 

rlomho et al. 

201)2. hereafter 


Papers I-IV). This cluster constitutes one of the best 
places to complete o ur objective, s ince i t harbours a large 
number of binaries: Ibevato et al.l 1 1991 ) found at least 5 
spectroscopic binaries and several radial velo city variables 
in Tr 16. The (very) preliminary results of iLevato et alJ 
Hffflllf . relying on poor quality and scarce data, were 
subsequently revised in the first papers of this series. 

In the present paper, we analyse an extensive set of 
spectra of one of the still poorly studied systems of Tr 16, 
HD 93161 (= Tr 16-176, CPD-SS^l) 1 . The compound 

1 The membership of this system to Tr 14, Tr 16 or Collinder 232 
has been subjeetto controversy throughout the literature. Here 
we follow IWalbornl Jl973l) in considering that it belongs to Tr 16. 
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Table 1 . Journal of the spectroscopic observations of HD 93161 A and B. Column 1 gives the instrument 
used and Column 2 lists the Heliocentric Julian Date (HJD) at mid-exposure. The next three columns 
present the phases (calculated from the Hell4471 orbital solution with e = 0) and the Hell4471 RVs of the 
primary and secondary components. The last three columns are similar, but refer to the averaged velocities 
(see Section 5). 


Instrument 

HJD—2 450 000 
(days) 

</>4471 

RVi 
(km s 

rv 2 

- 1 ) 

0avg 

RVi 
(km s 

RV 2 

- 1 ) 

HD 93161A 

B&C 

534.752 

535.758“ 

0.085 

90.0 

-105.4 





536.748 

0.318 

111.8 

-167.0 





537.741 

0.434 

81.8 

-83.9 





538.733 

0.550 

-20.1 

-20.1 





539.661 

0.658 

-151.0 

100.9 




BME 

1328.625 

0.760 

-178.5 

194.6 

0.761 

-156.9 

173.5 


1329.556 

0.869 

-112.2 

127.7 

0.870 

-117.3 

125.8 


1330.590 

0.990 

-6.7 

-6.7 

0.991 

-5.5 

-4.4 


1331.674 

0.116 

81.1 

-141.9 

0.117 

91.3 

-130.2 


1332.560 

0.220 

149.6 

-195.4 

0.221 

153.3 

-192.6 

FEROS 

1299.554 

0.367 

102.8 

-144.2 

0.368 

111.1 

-150.4 


1300.563 

0.484 

-4.0 

-4.0 

0.485 

-5.9 

-5.0 


1301.670 

0.614 

-100.9 

131.6 

0.615 

-102.1 

130.8 


1302.662 

0.729 

-157.6 

204.6 

0.731 

-150.3 

198.2 


1304.560 

0.951 

-54.7 

66.6 

0.952 

-44.7 

62.5 


1327.531 

0.633 

-98.9 

155.1 

0.634 

-102.6 

151.9 


1669.577 

0.562 

-55.6 

81.9 

0.563 

-58.4 

79.8 


1670.539 

0.675 

-139.7 

184.2 

0.675 

-143.3 

184.2 


1671.535 

0.791 

-151.8 

195.7 

0.791 

-149.2 

192.4 


1672.527 

0.907 

-92.7 

121.0 

0.907 

-95.7 

109.7 


1673.526 

0.023 

-3.6 

-3.6 

0.024 

-2.0 

-1.1 


2037.627 

0.528 

-70.7 

19.8 

0.528 

-7.9 

18.6 


2038.556 

0.636 

-126.1 

107.5 

0.636 

-120.5 

112.4 


2039.628 

0.761 

-159.8 

211.4 

0.761 

-156.5 

205.7 


2040.636 

0.879 

-79.0 

133.4 

0.879 

-93.0 

133.5 


2335.615 

0.314 

112.6 

-192.8 

0.314 

118.4 

-191.0 


2336.617 

0.431 

24.7 

-102.4 

0.431 

52.7 

-72.6 


2337.606 

0.547 

-70.3 

80.9 

0.546 

-64.7 

49.5 


2338.604 

0.663 

-111.7 

173.4 

0.663 

-114.5 

167.5 


2381.573 

0.679 

-138.5 

185.4 

0.679 

-139.4 

182.4 


2382.552 

0.794 

-156.8 

195.8 

0.793 

-154.3 

190.8 


2383.556 

0.911 

-75.7 

105.2 

0.910 

-79.1 

103.4 


2782.527 

0.486 

-2.2 

-2.2 

0.485 

-2.4 

-0.9 


2783.516 

0.601 

-84.2 

121.5 

0.600 

-87.9 

121.0 


2784.511 

0.717 

-157.4 

205.3 

0.716 

-154.4 

200.8 


3130.515 

0.109 

101.2 

-118.2 

0.108 

100.7 

-122.4 


3131.495 

0.224 

161.2 

-192.8 

0.222 

160.8 

-193.3 


3132.494 

0.340 

141.9 

-170.4 

0.339 

141.9 

-170.8 


3133.575 

0.466 

39.9 

-55.3 

0.465 

39.7 

-50.0 


3134.513 

0.576 

-65.4 

97.7 

0.574 

-68.1 

94.4 


3135.544 

0.696 

-146.9 

195.2 

0.695 

-147.1 

193.1 

HD 93161B 

BME 

1330.640 


-4.5 





FEROS 

2382.563 


-57.8 






3132.506 


-14.0 






a. No measurement is available for Hell4471, but the Hell 14686 line shows RVs of 109.2 and —161.2 kms 
for the primary and secondary, respectively. 
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formed by HD 93161 and HD 93160 actually harbours five 
stellar components, noted A -E, some o f which were found 
by speckle observations llMason et a 1.1 [Toolt l: HD 93160 
corresponds to the C component, located at 13” from 
the AB pair; component D is at 3” from HD 93160 and 
component E is at 8” of the AB pair. The latter two 
stars, separated by only 2", have been known for a long 
time (see e.g. |Walborr |H173|) and are generally r e ferred 
to as HD931 61~ TPennv iTotifil and IHowarth et all ill D97I') 
reported radial velocity (RV) variations in the UV spectrum 
of HD 93161. Unfortunately, the small separation between 
A and B, together with the low resolution of IUE, have 
not permitted to attrib ute these changes co nfidently to one 
member of the pair. IHowarth et alJ 1199 71 even suggested 
that b oth stars cou ld be short-period binaries. Previously, 
iLevato et al.1 1 1991 1 had provided an SB1 orbital solution 
for HD 93161, without specifying to which component this 
solution applied. A thorough study of the system was thus 
eagerly awaited. 

This paper is organized as follows: the observations 
and the data reduction are presented in Sect. 2 and the 
spectral features of HD 93161 A and B are discussed in 
Sections 3 and 4. In Sect. 5, we derive the first full orbital 
solution for HD93161A, whereas Sect. 6 deals with the 
analysis of its X-ray emission. As a next step, we examine 
the evolutionary status of the stars and we finally conclude 
in Sect. 8. 


2 OBSERVATIONS AND DATA REDUCTION 
2.1 Optical spectroscopy 

HD 93161 lies in a crowded region of the Carina Nebula, 
with HD 93160 just 13” west of the star. As noted in the 
introduction, HD 93161 itself was found to be a visual dou¬ 
ble system, composed of two stars with similar brightn esses 
separated by only 2”. As in IVreux fc Andrillatl il974f) , we 
define component A as being the one closest to HD 93160. 
With such a small separation between A and B, it is difficult 
to study these stars individually. During some observing 
nights, a poor seeing and/or the low angular resolution of 
the pointing camera did not enable us to clearly disentangle 
the close pair AB. Among our 45 spectra, only 3 can 
securely be attributed to the sole B component. The main 
contributor to the remaining 42 spectra is intended to 
be the A component but these observations are usually 
contaminated, to a variable degree, by the B companion. 
This is marked in the spectra by the presence of a third, 
weak component. 

We have observed HD 93161 over a period of 7 years 
with various instruments. A first set of 6 medium resolution 
spectra covering the wavelength range 3850 — 4800 A was 
gathered in 1997 with the ESO 1.5 m telescope equipped 
with a Boiler & Chivens (hereafter B&C) Cassegrain 
spectrograph. The data were obtained with a holographic 
grating (2400 lines/mm, ESO grating #32) providing a 
reciprocal dispersion of 32.6 A/mm. The detector was a 
thinned, UV flooded Loral-Lesser CCD (ESO #39). The 


slit width was set to 220 /im corresponding to 2” on the sky. 
The spectral resolution as measured on the helium-argon 
calibration spectra is 1.2 A, corresponding to a resolving 
power of 3600. Typical exposure times were of the order 
of 10 minutes and the average signal-to-noise ratio was 
about 250. Part of the BfeC spectra were affected by a 
fringing pattern llTuratto et alJ Ill997 i that occurred over 
the wavelength range 4050 — 4250 A in this instrument con¬ 
figuration. Given the variability of the fringing pattern and 
in order to avoid amplification of the fringes in the stellar 
spectra, the data were not flat-fielded. All the reductions 
were performed using the MIDAS software developed at 
ESO. Note that only B&C spectra were sky subtracted. 

Six high resolution spectra of HD 93161 were taken in 
May-June 1999 with the Bench-Mounted Echelle spectro¬ 
graph (BME) fed by the 1.5 m Ritchey-Chretien telescope at 
CTIO. These data covered the spectral range 3750 — 5800 A. 
The typical S/N ratio was ~60 at 5430 A for exposure times 
of one hour. The detector used was a Tek CCD with 2048 
x 2048 pixels of 24 /im x 24 /im and the resolving power 
was 45 000. The BME data were reduced using the IRAF 2 
package, following the recommendations of the BME User’s 
Manual. A first rectification of the extracted orders was 
carried out with the projector flat exposures. The spectra 
were then normalized by fitting a low-order polynomial to 
the continuum. 

Between 1999 and 2004, we obtained 33 additional 
high resolution spectra at the La Silla Observatory with the 
Fiber-fed Ex tended Range Optical Spectrograph (FEROS, 
iKaufer et all 1999). This spectrograph was installed un¬ 
til October 2002 at the ESO 1.5 m telescope, and was 
subsequently moved to the 2.2 m telescope. The FEROS 
spectra cover the wavelength range 3750 — 9000 A with 
a resolving power of 48 000. The detector was a 2k x 4k 
EEV CCD with pixel size 15 (im x 15 /im. The typical 
exposure time was 10-15 minutes and the S/N ratio at 
5650 A is about 125-225. The FEROS spectra were reduced 
using the appropriate context of the MIDAS e nvironment, 
toget her with an improved reduction pipeline ISa.na et a.lJ 
1200311 . They were normalized by fitting the continuum with 
a polynomial of degree ranging from 2 to 6. Most of the 
unavoidable fringes affecting the red part of the spectra 
were corrected by flat-fielding, but there are a few residuals, 
e.g. near 6640 A. 


2.2 X-ray data 

The core of the Carina Nebula wa s observe d sev eral times 
by the XMM-lVeuiton observatory d.Tansen et a.lJ 1200I V We 
will focus here on the five datasets acquired in 2000 and 
2001. The first two datasets, centered on ?; Carinae, were 
obtained in July 2000 during satellite revolutions #115 and 
#116. The three other observations, centered on WR25, 
were taken one year later during Revs. #283, #284, and 


2 IRAF is distributed by the National Optical Astronomy Obser¬ 
vatories. 
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Table 2. Journal of the XMM -Newton observations of HD 93161. 
The first column gives the revolution during which the observa¬ 
tions were taken, Column 2 the Julian Date (JD) at mid-exposure, 
Column 3 the phases (calculated from the average He I orbital so¬ 
lution, see Sect. 5), and the last three columns the exposure times 
inks for each instrument. For more details, see iRaassen et alJ 

boofl . 


Rev. 

JD—2 450 000 
(days) 

4>awg 

Exp. Time (ks) 

MOS1 MOS2 PN 

115 

1751.92 

0.18 

33.8 

30.6 

31.7 

116 

1753.56 

0.37 

11.2 

8.3 

9.4 

283 

2086.00 

0.17 

36.7 

36.7 

34.6 

284 

2089.05 

0.53 

42.1 

42.1 

39.6 

285 

2090.91 

0.75 

37.4 

37.4 

34.9 


#285. A preliminary stud y us ing only th e first two obse rva- 
tions was presented bv lAlbacete Colombo et alJ (2003). A 
thorough analysis of the field relying on t he e ntire observa - 
tional dataset has been performed hv lAntokhm et ailEopa. 

HD 93161 was observed serendipitously during all five 
XMM-Newton pointings (see Table 0 The source was 
recorded with the three European photon imaging Camera 
(ep ic) detectors (PN. MOSl and M OS2; for more details, 
see lStriider et af]l200lt [Turner et alJl200ll) . The data were 
processed with the XMM-lVeiofon Science Analysis Software 
(SAS, version 5.4.1). A check on proton flares resulted in 
the deletion of part of the exposure time. More details 
on the pip elin e proc essing of these data can be found in 
iRaassen et al.l ( 20oJ l 3 . 

3 THE SPECTRUM OF HD 93161A 

Up to now, the visible spectrum of HD 93161 was essentially 
stu died through photographic observations. Following 
I Walborn] ; • 1973 both stars of HD 93161 have identical 
spectra and apparent magnitudes. Walborn further at¬ 
tributes an Q6 .5V((f)) type to e a ch sta r. The following 
year, however, IVreux fc Andrillatl f 1974 1 discovered that 
HD93161A was cooler than HD93161B. Unfortunately, 
no formal spectra l type was assigned to HD 93161A by 
IVreux fc Andrillatl lll974lf . 

We have investigated the spectrum of HD 93161A 
through high resolution echelle spectroscopy. Such a spec¬ 
trum is shown in Fig.0 the lines are clearly double, showing 
that this star is indeed a spectroscopic binary. Beside the 
lines of the Balmer series, the spectrum reveals many 
absorption lines of Hel, Hell, Nill, Silll, Siiv, Civ and 0 
that are typical of O stars. It also reveals a number of diffuse 
interstellar bands (DIBs), as well as interstellar absorptions 
due to He113889, Nai, Can, CH and CH + . There are only 
few emission lines in the spectrum of HD 93161A. Most of 
them are sharp nebular lines, associated with the Carina 
Nebula. We may note that the most massive component 

3 In R.aassen et alJ 1200 . 4 ) , note that Column 4 of Table 5 actually 
gives the Julian dates in the format JD—2 450 000, not the MJDs. 


exhibits Cm 15696 in emission 4 . Throughout this paper, 
we will refer to the more massive star as the primary, or 
HD 93161A1. 

Only the least blended and higher quality spec- 
tra are suitable for the spectral classificatio n. F ollowing 
IConti fc_Alschuie3lll97lh : ICont,i fc Frostldl977ll and lMathvsl 
ll 984 ll 9891) . we adopt, the log [MF (4471)/EW (4542)] clas¬ 
sification criterion to determine the spectral types of 
the stars. For this purpose, we fitted these lines with a 
combination of several gaussians. This procedure enables 
to measure the radial velocities (RVs) rather securely, but 
the presence of small additional components (due to the 
contamination by the neighbouring HD 93161B) renders the 
determination of accurate equivalent widths (EWs) rather 
awkward. The measured ratios are thus rather dispersed, 
resulting in an approximate 07-7.5 type for the primary 
and an 07-8.5 type for the secondary at 1-a dispersion. To 
refine this classification further, we h ave compared their 
spectr um to the OB star atlas of I Walborn fc Fitzpatrickl 
1 1996 ). For the primary star, Hell4471 (resp. 4026) appears 
only slightly stronger than He ill 4542 (resp. 4200). The 
reverse situation occurs for the He ill 4143, 4388 lines, 
which are weaker than the He nil 4200, 4542 lines. The 
spectrum o f H D 93161A 1 is t hus similar to that of HD 48279 
dWalbom fc Fitzpatricklll99oi) . indicating an 08 type for 
this star. On the other hand, the secondary seems cooler: 
He 114471 (resp. 4026) is much stronger than He ill 4542 
(resp. 4200), while Hel 14143 is weaker than Hel 14388 
but has a strength similar to those of He nil 4200, 4542. 
HD93161A2 thus resembles HD 46202, an 09 star. These 
spectral types are slightly later than those derived from 
Conti’s criterion. This difference is probably due to the 
contamination by the B component but we can not totally 
exclude that the stars might be one subtype earlier. 

For the luminosity classification, we rely on the 
log [EW( SiIV A 4089)/EW(He IA 4 143)1 criterion intro¬ 
duced bv IConti fc Alschulerl dl97lll . For the same reason as 
before, the measured values of this ratio are rather scattered 
but they clearly indicate a main sequence classification. 
Note that the log [EW (HeIIA 4686)1 +log \EW (Rel A 4388)1 
luminosity criterion introduced by iMathvsl ( T988I . Il989h 
was not used here: it relies on absolute EWs, which are 
difficult to determine because of the variable contamination 
by HD 93161B. The main sequence classification is further 
supported by the comparison of our spectra with the atlas 
of standards since, for both stars, Si IV1 4089 is quite faint 
whereas the He ill 4686 line appears rather strongly in 
absorption. We thus conclude that HD93161A is most 
probably an 08V+09V binary. 


4 The fact that only this star exhibits the C III 15696 line in emis¬ 
sion helps for the identification of the primary star. In addition, 
the presence of a second emission peak in this line clearly indi¬ 
cates pollution by HD93161B, which also displays this line in its 
spectrum. It thus enables us to estimate the degree of contami¬ 
nation by HD 93161B in each spectrum of HD 93161A. 
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Figure 1. Spectrum of HD93161A as observed with the FEROS spectrograph on HJD 2453 135.544 (0 ~ 0.7). At this orbital phase, 
the spectrum of the most massive star is shifted towards the blue. The main lines are indicated at their rest wavelengths. 


4 THE SPECTRUM OF HD 93161B 

A FEROS spectrum of HD 93161B is shown in Fig. [31 This 
spectrum is similar, but not identical to that of HD 93161A. 
HD93161B actually exhibits a few more emission lines, 
e.g. the weak S IV 14485, N III 114634, 46 41 lines. As already 
remarked by lYreux fc Andrillatl ill 9741) . the star appears 
slightly hotter than its A neighbours, since the Hell4471 
line is weaker than the He ill 4542 line. 

Among the available spectra of HD 93161B, 
we again rely only on the high-quality FEROS 
data for the spectral classification. We used the 
same classification criteria as above: we deter¬ 


mined log [£W(HeI A4471) / £W(HeII A4542)] 

—0.14, suggesting an 06.5 spectral type, while 
log [EW (Si IV A 4089) / EW (He IA 4143)] < 0.1, corre¬ 

sponding to a main sequence classification. We have also 
compared the spectrum of HD 9316 1B with those of the 
spect ral atlas of OB standards dWalborn fc Fitzpatrick! 
Il990l) . In the spectrum of HD93161B, the He nil 4200, 
4542 lines are much stronger than the He ill 4143, 4388 
lines whereas Hell 4026 is stronger than He ill 4200. This 
confirms the 06.5 type for this star. The relative weakness 
of the Si IV 14089 line indicates a main sequence classihca- 
tion, again in agreement with Conti’s criterion (which in 
principle does not apply to stars hotter than 07). Moreover, 
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A (A) 

Figure 2. Spectrum of HD 93161B as observed with the FEROS spectrograph on HJD 2 453 132.506. The same lines as in Fig. ITI are 
identified. 


the N ill 114634, 4641 lines exhibit only a weak emission 
while the He ill4686 line is in absorption (with a lower 
intensity than that of the He ill 4542 line) , as is ty pically 
seen in some main sequence stars llWalborn fc Fitzpatrick! 
il990t ). The presence of these latter lines require the addition 
of an ((f)) tag to the classification: HD93161B is thus an 
Q6 .5V((f ) ) sta r, in agreement with the classification of 
IWalbornl 1 1973 ). 

The bottom part of Table 0 gives the measured RVs 
of the He 1 14471 line on our 3 spectra of HD93161B. 
While we find no obvious evidence for a companion in the 
spectrum, the observed RVs of all lines of HD 93161B are 


clearly variable. Except perhaps for the lower quality BME 
data, this is most likely not due to contamination by the A 
neighbour. Indeed, the contemporaneous spectra of the A 
star show two widely separated components, which are not 
(even slightly) detected in the B spectrum and could thus 
not hamper the RV determination. These RV variations 
might actually be caused by orbital motion in a binary. 
Unfortunately, the small number of B spectra does not 
enable us to further analyse the variability of this star, 
and especially to determine if these variations are periodic. 
Additional uncontaminated, high S/N data of HD93161B 
would certainly allow to investigate the multiplicity of this 
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Table 3. Circular orbital solutions for HD93161A as deduced from different lines. K\ and K 2 yield the amplitude of the RV curves 
of the primary and secondary, respectively, whilst 71 and 72 correspond to the apparent systemic velocities of the two components. To 
account for a possible velocity gradient in the line formation region, we allowed 71 and 72 to be different. Quoted uncertainties are the 
1-er error bars. s y /s x gives the ratio of the secondary to primary uncertainties on the RVs. 


Lines 

s y/ s x 

p 

(days) 

K! 

(km s — *) 

k 2 

(kms -1 ) 

71 

(kms -1 ) 

72 

(kms -1 ) 

rms 

(km s -1 ) 

He 1 14026 

0.8 

8.5665 

149.7±3.4 

201.24=4.6 

5.14=4.0 

13.04=5.0 

14.8 

He 1 14471 

1.0 

8.5662 

151.54=2.4 

202.24=3.2 

-1.44=2.5 

2.24=3.1 

12.4 

He 11 14542 

1.4 

8.5664 

150.8±4.8 

198.14=6.4 

4.7±7.7 

2.94=9.6 

20.1 

He 11 14686 

1.1 

8.5663 

150.84=2.4 

201.14=3.2 

9.14=2.5 

3.04=3.0 

12.4 

He 1 14921 

0.8 

8.5662 

148.94=2.3 

202.34=3.1 

-1.14=2.1 

1.54=2.4 

9.4 

He 1 15016 

0.8 

8.5662 

155.34=1.6 

201.0±2.1 

0.64=1.3 

-0.24=1.4 

6.6 

Hell 15412 

1.3 

8.5666 

150.84=4.1 

192.94=5.2 

9.0±7.5 

-0.24=9.5 

20.8 

015592 

1.2 

8.5664 

155.24=2.3 

201.34=3.0 

0.74=2.3 

-4.24=2.7 

10.4 

Civ 15801 

1.2 

8.5663 

153.04=2.7 

198.34=3.5 

—2.54=2.9 

-5.44=3.5 

11.3 

Civ 15812 

1.4 

8.5663 

155.14=2.4 

197.94=3.1 

-5.14=2.5 

-2.84=3.0 

10.0 

He 1 15876 

0.9 

8.5663 

152.54=2.4 

202.24=3.2 

3.94=2.4 

3.14=2.8 

9.8 

He 1 17065 

0.9 

8.5663 

155.84=2.2 

201.44=2.9 

2.94=2.0 

4.14=2.2 

8.4 


star. 


5 ORBITAL SOLUTION FOR HD 93161A 

As mentioned before, we have determined the RVs of the 
lines by fitting them with a sum of gaussian profiles. In order 
to get a reliable orbital solution, we considered 12 different 
lines: He ill 4026, 4471, 4921, 5016, 5876, 7065; He nil 4542, 
4686, 5412; 015592 and Civ 115801, 5812. We present in 
Table Q the RVs measured for the Hell4471 line. Note 
that we have also measured the RVs of the most prominent 
interstellar lines: for our FEROS observations, the average 
RV of the CH + 14232 line is 5.3 ± 0.4 kms -1 , whereas 
for the B&C data, the mean RV of the Ca ill3933 line is 
—35 ± 15 kins' 1 and for the BME data, the mean RV of the 
red component of the Ca 11 13933 line is —0.1 ± 2.9 kms -1 . 
We thus emphasize the very good temporal stability of 
the FEROS spectrograph, whose spectra constitute the 
majority of our sample. 

Up to now, two orbital solut ions for HD 93161 have 
been published in the literature. iLevato et ali j 1994.) have 
calculated an SB1 solution for HD 93161, but without 
specifying t o which component (either A or B) it refers. 
iLuna et alj j2003jl_ failed to pr ovide any improvement to the 
results of lLevato et alJ (4 99 if) 5 . If these authors measured 
the RVs on combined spectra of HD 93161, as can be 
guessed from their description, it seems rather strange that 
they detected only one component (thus deriving an SB1 
solution) since the presence of multiple components is easily 
spotted even in our medium resolution B&C spectra. In the 

5 Note that in lLuna et al 1 <l2003l) , the figure showing the ‘new so¬ 
lution’ for HD 93161 is most probably erroneous. Moreover, in the 
same bulletin, the authors also provide an SB1 solution for Tr 16- 
104 and claim to detect no RV variations for Tr 16-110, although 
both sta rs were already known as triple systems prior to th is pub¬ 
lication jRauw et aljl200ll; lAlbacete Colombo et al. 1 12002'). This 
may legitimately cast some doubt on the quality of their data. 


Table 4. Orbital and physical parameters of HD 93161A as de¬ 
duced from the mean RVs for a circular orbit. The usual notations 
are used, Rrl being the Roche lobe radii of the stars and To cor¬ 
responding to the conjunction with the primary star in front. The 
error on the period is calculated from the observational time base 
(see §5.1) and all other errors are estimated 1 -cr deviations calcu¬ 
lated for a fixed s y /s x ratio and a fixed period (the ones that give 
the lowest residuals). Note that the apparent systemic velocities 
(71 and 72 ) of each line were subtracted prior to averaging. 


P 

8.56634=0.0040 days 

T 0 (HJD 

3001.098 

-2 450 000) 

4=0.011 

Sy / S x 

0.9 

mi / m 2 

1.31±0.02 


152.94=1.7 kms -1 

k 2 

200.64=2.2 kms -1 

71 (after sub.) 

0.024=1.4 kms -1 

72 (after sub.) 

0.054=1.6 kms -1 

ai sin i 

25.94=0.2 R 0 

CL2 sin i 

33.94=0.4 Rq 

mi sin 3 i 

22.24=0.6 M 0 

m2 sin 3 i 

17.04=0.4 Mq 

Rl RL 

0.40264=0.0015 

a l + a 2 


r %t. 

0.35574=0.0014 

ai+a -2 


rms 

7.7 kms -1 


rest of this paper, we will not discuss these analyses further. 


5.1 Period determination 

Once the RVs are measured, the next step is to determine 
the period of the orbita l motion. To this aim, we have 
applied the method of lLafler &: Kinmanl l) 19651 here- 
after LK) as we ll as a Fourier-type analysis (method of 
iHeck hereafter HMM, see also the comments 

hv lOosset et a.1.112001 If . Since confusion could be a possible 
issue, these period search algorithms were first applied to 
the time series of the absolute RV differences \RVi — RV 2 |: 
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Figure 3. Radial velocity curves corresponding to the averaged 
RVs of HD 93161A and assuming e = 0.0 (see TableTH. The open 
and filled circles correspond to the observed mean RVs of the pri¬ 
mary and secondary, respectively. Note that we used a reference 
frame with zero systemic velocity (see text). 

for each of the lines quoted above, a period of about 4.3 days 
was then detected. Using this preliminary period and the 
RVs of the C ill 15696 line, we then identified the primary 
and secondary, and applied the period search on the actual 
differences RV 1 — RV 2 . Both the Fourier analysis and the 
LI< algorithm agree on a value of P = 8.565 days for all lines. 

Our dataset consists of 42 spectra spread over 2600 days 
for the He 1 14471 and He 11 14686 lines (where there exist 
B&C data in addition to the FEROS and BME obser¬ 
vations) and of 36 spectra spread over 1836 days for the 
other lines. These observational time bases correspond 
respectively to natural peak widths Au of 0.00038 d~ x 
and 0.00054 d -1 . The full width at half maximum of the 
actual peaks in the periodogram is in agreement with these 
theoretical values. We adopt an uncertainty of one tenth of 
the peak’s width, i.e. P = 8.565 ± 0.003 or 0.004 days, the 
uncertainty depending on whether or not the B&C data are 
considered. 


5.2 Orbital elements 

To de rive t h e orb ital elements, we used a modified version 
of t he 1 Wolfe et al. (i 1 1 9671) method. The algorithm, presented 
by ISane^t alJ l2003ih allows a weighting that differs for 
the secondary and primary RVs. In the present paper, 
we have adopted a period value and a relative (primary 
vs secondary) weight that yield the lower residuals. In 
addition, we accounted for the quality of the measurements 
by attributing lower weights to the observations where 
the components are blended (i.e. only one measurement is 
available for both stars), noisy (e.g. BME and B&C data), 
or heavily contaminated by HD93161B. 


At first, we fitted the RVs from each line with a 
general, eccentric solution. However, the computed ec¬ 
centricity e turned out to be very small (< 0.03), with a 
null eccentricity within less than two sigma. Therefore, a 
circular orbit was assumed for HD 93161A. The results of 
these non-eccentric fits are displayed in Table [3 for each 
line. They are compatible with each other, within the errors. 

Finally, we tried to get more accurate results by 
averaging the RVs from different lines. This was done after 
subtracting the apparent systemic velocities derived for each 
line from the individual fits. Note that the lines yielding 
the largest residuals (i.e. He nil 4542, 5412 and Hell 4026) 
were not considered. The lower resolution B&C data were 
also discarded in order to have a more homogeneous set of 
RVs. Averaging all other observations allows to obtain an 
rms residual of 8 kins' 1 . This best solution, including a 
reappraisal of the period (Pfi na i = 8.5663 ± 0.0040 days), 
is presented in Table 21 and is shown in Fig. [3 it will be 
adopted throughout the rest of the paper. 

The minimum masses m sin 3 i derived for the primary 
and secondary are 22.2 M@ and 17.0 Mq, respectively (see 
Table 2J • Comparing these values to typical masses of 08V 
and 09V stars allows to estimate the inclination of the 
system. However, one must be careful when doing so, since 
the analysis of the binaries in the Tr 16 cluster (HD 93205, 
Paper II and Tr 16-104, Paper III) has shown that the 
typical masses pre s ented by well-known references like e.g. 
iHowarth & Prinial lll989l) are generally too large compared 
to the masses derived from the study of binary systems. 
We therefore decide to ad opt the v alues determined for 
detached eclipsing systems dCiesll2003l) : ~22 M 0 for 08V 
stars and ~18 M@ for 09V stars. Such values are very 
close to the minimum values found for HD 93161A 1 and 2, 
indicating a large inclination, most probably > 75°. With 
such an inclination, the system should undergo eclipses. 
Unfortunately, no photometric study of HD 93161 has been 
conducted up to now: additional observations are thus 
needed. They will help to constrain further the physical 
parameters of HD 93161A. 


6 XMM-lVeioton OBSERVATIONS OF HD 93161 
6.1 Spectral analysis 

Due to the proximity of HD 93160, the spectrum of 
HD 93161 was extracted from a circular region with a 
radius limited to 10”, as is shown in Fig. 21 Due to a 
poor calibration of the Point Spread Function of the epic 
instruments (especially PN) at large off-axis angles, we 
extracted MOS and PN spectra in the extended source 
mode, i.e. without any EEF (encircled energy fraction) 
correction. For the chosen region, we estimate that this 
correction amounts to a factor ~ 1.8. The background 
spectrum was obtained by selecting events from a circular 
region spatially offset from the source (see Fig. 21- 

The X-ray data were then analysed with the xspec 
software (version 11.2.0). For each observation, the epic 
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Table 5. Best-fitting models and X-ray fluxes at Earth for each XMM -Newton observation of HD 93161. 
The fitted mo del has the for m wabsCTV^J.)*(mekal(kT\ )+wabs(7V H )*mekal(kT2)), with wabs(iV^. )= 4.5 X 
10 21 cm 2 ID ini as ± Sinavfll994l ). Quoted fluxes are in the 0.4—10. keV energy range and are not corrected 
for EEF (see Sect. 2.2). The unabsorbed fluxes /)(. n ' ll>s are dereddened only considering the interstellar 
absorbing column. The quoted errors corresponds to 1-cr deviations. 


Rev. 

^avg 

kTi 

keV 

Nf 

10 22 cm -2 

kT 2 

keV 

xS(dof) 

f abs 

J X 

(10 -13 

funabs 

J X 

erg cm -2 s _1 ) 

115 

0.17 

0.26 ±0.02 

0.43 ±0.13 

0.71 ±0.11 

1.14 (132) 

0.85 

4.91 

116 

0.37 

0.29 ±0.03 

0.44 ± 0.66 

0.93 ±0.67 

1.26 (37) 

0.77 

4.46 

283 

0.17 

0.27 ±0.03 

0.45 ±0.14 

0.67 ±0.08 

0.87 (106) 

0.92 

4.67 

284 

0.53 

0.29 ±0.02 

0.44 ±0.13 

0.71 ±0.12 

0.89 (122) 

0.96 

4.72 

285 

0.73 

0.29 ±0.01 

0.86 ±0.16 

0.79 ±0.11 

1.18 (139) 

0.89 

4.57 


HD93161 : all EPIC in 0.4-12 keV 



Right Ascension (J2000) 

Figure 4. Combined EPIC image showing HD 93161, HD 93160 
and their neighbourhood. The bright X-ray source to the top right 
is the core of the Tr 14 cluster, that contains notably Tr 14-1 and 
HD 93129AB. The source and background regions used to extract 
the X-ray spectrum of HD 93161 are also shown. Note that this 
image was exposure-corrected. 


spectra were fitted simultaneo usly by a two-temperature 
therm al plasma mekal model (IMewe et a.l .111 9851 : IK austral 
1992) of the form wabsi nt *(mekali+wabs*mekal 2 ) where 
wabs int was fixed to th e int erstellar absorption determined 
by iDiplas fc Savagd ll 19941 N^ t = 4.5 x 10 21 cm -2 ). 
The second wabs model is included to account for the 
possibility of additional absorption for the hard component 
of the spectrum. The abundances were set to solar. The 
best-fitting model parameters are presented in Table 0 
and the spectra from Rev. #284 are shown in Fig. 0 along 
with their best-fitting model. The spectral parameters and 
the fluxes are fairly stable throughout the five datasets. 
The parameters are quite typical for a single O-type star, 
with kTi ~ 0.3 keV along with a slightly higher temper¬ 
ature component. Our resul ts are in agree ment w ith the 
preliminary conclusions of lAlbacete Colombo et’afl (2003) 
concerning the first two datasets. 



o.e i 2 s 

channel energy (keV) 


Figure 5. EPIC spectra of HD 93161 obtained for Rev. #284. In 
the top panel, the upper and lower data correspond to the PN and 
MOS spectra, respectively, and the solid lines yield the best-fitting 
absorbed two-temperature model (see Tabled. The bottom panel 
shows the contributions of individual energy bins to the x 2 of 
the fit. The contributions are carried over with the sign of the 
deviation (in the sense data minus model). 


6.2 Variability in the X-ray domain 

We also investigated the variations of the X-ray emission 
of HD 93161, first on a short time scale (within each point¬ 
ing) and then on the long-term range (between the differ¬ 
ent pointings). For each dataset, the source and background 
events were extracted using the same regions as those de¬ 
fined for the spectrum analysis. After taking into account the 
flaring intervals (see Sect. 2.2), these event lists were used to 
construct background-subt racted lightcurve s that were anal¬ 
ysed by various methods llSana et alJ 12004 ). No significant 
short-term variability was consistently detected in the three 
instruments. 

As a next step, we searched for variations between the 
different pointings. To this aim, we compared the mean 
count rates as well as the fluxes determined by spectral fit¬ 
ting (see Tabled: no significant long-term changes from one 
pointing to another were detected. 


6.3 Wind interactions in HD 93161? 

We found no clear evidence that the components of 
HD 93161 are interacting in any manner. The visible spec¬ 
trum presents only few emission lines and, in particular, 
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no Ha em ission is seen that might arise in a wind colli¬ 
sion zone llThallerlll997t ISana et all 12003) . Moreover, the 
minimum Roche lobe radii, Rrl sini, determined for the 
primary and secondary of HD93161A are 24.1±0.2 and 
21.3±0.2, respectively (see Table 0. These large values con- 
trast with the small typical radii (< 10 R(t 0 of Q8-9V stars 
Jlfow^di ^Prima[_^89||_|Morrell et, alj 1200 it iRauw et alJ 
l2QtHnGie^l200^1SaimetaLll20o'5lJl . indicating that this ob¬ 
ject is a well detached binary system without any ongoing 
Roche lobe overflow. In addition, HD 93161 is not a par¬ 
ticularly bright X-ray source, since its unabsorbed luminos¬ 
ity amounts only to ~ 6.3 x 10 32 erg s _1 (after EEF cor¬ 
rection and for a distance of 2.5 kpc, see next section), or 
log(Lx/Tboi) = —6.6 (using the Lbol determined in the next 
section). 

We can compare this observed X-ray luminosity to the 
emission expected for a wind collision in HD 93161A. First, 
we determine the mass-lo ss rates M a nd ter minal velocities 
Voc using the recipes of lYink et alJ ll200ft 1200 if) and the 
probable radii, masses, temperatures, and bolometric 
luminosities of the stars (see Sects. 5 .2 and 7). Using the 
formalism of IPittard fc Stevensl ll2002l) . we then find that, 
in the case of HD 93161A, a radiative wind collision would 
lead to an X-ray emission with Lx ~ 2 x 10 34 erg s _1 . 
However, the physical parameters of HD 93161A (strengths 
of the winds, separation between the stars) should rather 
lead to a slightly adiabatic shock: the X-ray emission 
associated with such a collision could thus be significantly 
lower than in the radiative case. In fact, the total luminosity 
derived after applying the ‘ canonical’ Lx-L^o i relation (see 
iBerghofer et alJll997l . and [Sana et al J l2005ah to the three 
massive stars of HD 93161 is ~ 3.4 x 10 32 erg s _1 , indicating 
that the observed X-ray emission does not display a large 
excess attributable to a wind interaction. This conclusion 
is further supported by the lack of hard X-ray photons 
and the constancy of the X-ray emission. We might note 
however that our X-ray observations do not sample the 
whol e orbi tal cyc le: as was the case for CPD—41°7742 
( Sana et al . 2005bl) , the signature of a wind collision could 
be limited to a few crucial phases. 

Finally, iBenaglia fc Koribalskil (' 20041 1 observed the 
Carina nebula with ATCA at 3 and 6 cm aiming at the 
detection of non-thermal radio emission but HD 93161 was 
not detected in this survey, contrary to what would be 
expected if there was a wind collision between HD 93161A 
and HD 93161B. Therefore, from the available data, we 
conclude that there is probably no strong wind interaction 
in HD 93161. 


7 FUNDAMENTAL PARAMETERS 

Using the temperature scale of lMartins et al.l 1200 21. we can 
convert the spectral types derived above into effective tem¬ 
peratures of 38200 ± 700 K, 35250 ± 700 K and 33700 ± 800 K 
for HD93161B, HD93161A1 and HD93161A2 respectively. 
The quoted uncertainties correspond to half a subclass for 
the components of HD93161A and to the range of temper¬ 
atures for the 06.5 V class according to Martins et al. 

In order to derive the luminosities, we need the ob¬ 


served magnitude, the reddening, the brightness ratios be¬ 
tween the different components, the bolometric correc ti ons 
and an esti mate of the distance. The WEBDA 6 iMermilliodl 
Il988[ll992l) database lists the results of several photometric 
studies that provide the apparent V magnitude of HD 93161 


JWalborn! 1973 l t Feinstein et al. 1972f: 

Kaltcheva &; Georeiev 

199,4 Forte 1978: Massey & Johnsonl 

19931). On average, we 


find V = 7.84 ± 0.02 for the entire HD 93161 system. On 
the other hand, lYazauez et alJ i l996i) quote V = 8.56 and 
V = 8.60 for HD 93161A and HD 93161B respectively. These 
magnitudes yield V = 7.83 for the entire system (compo¬ 
nents A + B ) in good agre ement with the mean value here- 
above. iMason et ahl ( 1998 1 found V = 9.0 for HD93161A 
and V = 9.1 for HD93161B. If these values were correct, the 
total magnitude of HD 93161 should be 8.3 in contradiction 
with the other observations. In the following, we therefore 
adopt apparent magnitudes of 8.56 ± 0.02 and 8.60 ± 0.02 
for HD 93161A and B respectively. 

The ( B — V) colour indices of HD 93161A and 
HD93161B were found to be 0.20 and 0.23 respectively 
iVazauez et al.l ll99fj). Comparing with typical colour in¬ 
dices of O stars lISchmidt-Kalen 1198211 . we hence adopt 
E(B - V) = 0.52 and E(B - V) = 0.55 for HD93161A 
and HD 93161B. 

To estimate the relative brightness of the components 
of HD 93161A in the visible spectral range, we evaluate the 
ratio of the EWs of the primary and secondary components 
for a number of lines in our 2004 FEROS spectra (the least 
contaminated by the B component). We use the He ill 4471, 
5876, He II14542 and 015592 lines and we compare the EW 
ratios to the values fo r typical EWs evalu a ted fr om th e com¬ 
pilatio n provided by IConti fc Alschuieil (Il97lll and IContil 
( 1974 1. In this way, we estimate a brightness ratio of 1.7±0.3 
between the primary and secondary of HD 93161A. This re¬ 
sults in apparent magnitudes of 9.06 ± 0.07 and 9.64 ± 0.12 
for the primary and secondary components respectively. 

The distances to the clusters Trumpler 14, Trumpler 16, 
Collinder 228 and Collinder 232 are still poorly known. How 
controversial the issue of the distances and even the actual 
existence of these clusters is, can b e illustrated by a brief 
overview of the literature. IWalbornl ( 1995 1 summarised the 
situation as follows: Tr 14 is younger than Tr 16, but both 
clusters lie at the same distance; neither Collinder 228 nor 
232 are genuine clusters, the former is part of Tr 16, whereas 
the latter contains stars of Tr 14 and Tr 16. The ambiguities 
on the distances of Tr 14 and Tr 16 are at least partly due to 
the uncertainties on the extinction towards these clusters. 
Using UBV photometry. iMassev fc Johnson Il99 3;l inferred 
a distance of 3.2 kpc for Tr 14 and Tr 16. Assuming a nor¬ 
mal exti nction l aw an d b ased on Str omgren and H /3 pho¬ 
tometry, iKaltcheva fc Oeorgievi (j 1993 1 derived distances of 
individual stars in Tr 16. These authors reached a radically 
different conclusion: since our line of sight towards Tr 16 is 
almost parallel to the molecular cloud ridge in the Carina 
spiral arm, the stars might not belong to an actual cluster 
but could instead lie along an extended star formation re- 
gion w hich is project ed on a small area on the sky. Recently, 
ICarraro et ali ll2004tl presented a study of Tr 14, Tr 16 and 
Collinder 232 based on UBVRI photometry. These authors 


available at the url http://obswww.unige.ch/webda/ 
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argued that a unique reddening law is not appropriate to 
study the whole region and they estimated different values 
of the selective extinction Rv = Av/E(B — V) for each clus¬ 
ter. By comparing their observations with empirical main- 
sequences, they derived distances of about 2.5 kpc for Tr 14 
and Collinder 232, but a much larger distance of about 4 kpc 
for Tr 16. The latter result is clearly at odds with previous 
distance determinations. In parti cular, the studies of two 
eclipsing binari es in Tr 1 6 (T r 16-l. lFrevhammer et alJl200lL 
and Tr 16-104, iRauw et al ] {200i allowed to obtain abso¬ 
lute lumino sities of t hese stars. A dopti ng bo lometric correc¬ 
tions from iHumphrevs fc Mac Elrovl ( 19841) and assuming 
Rv = 3.3, a distance of about 2.5 kpc is inferred for these bi- 
na ries, i n exc elle nt ag reement with the Tr 16 distance found 
bv iTapia et al.l ( 2003 '). In the following, we shall therefore 
adopt a distance of 2.5 kpc for HD 931 61. 

For Collinder 232, ICarraro et alJ (l2004h infer Ry = 
3.73 ± 0.03, whereas they estimate Rv = 3.83 ± 0.33 for 
HD 93161 and Rv = 3.60 ± 0.27 for its close neighbour 
on the sky HD 93160. Adopting an average value of Rv = 
3.72 ± 0.16, we obtain absolute magnitudes of —4.86 ± 0.11, 
-4.29 ± 0.15 and -5.44 ± 0.09 for HD93161A1, A2 and 
B res pectiv el y. F inally , adop ting the bolometric corrections 
from [Humphreys fc Mac Elrovt 1198 40, we obtain bolomet¬ 
ric luminosities of logl/j^/Lg = 5.25 ±0.04, log L^/Lq = 
4.93 ± 0.06 and logL^ ol /L 0 = 5.53 ± 0.04. 

Using these values, we have plotted in Fig.[S] the compo¬ 
nents of HD 93161 in the Hertzspru ng-Russe ll diagram along 
with the evolutionary tracks from ISchaller et all ( 49921 ') for 
solar metallicity. This diagram suggests that the stars have 
an age of 3-4 Myr and are thus slightly evolved off the zero 
age main sequence (ZAMS), unlike other stars of the Tr 16 
cluster, e.g. the components of Tr 16-104 whi ch w ere found 
to lie very close to the ZAMS feauw et al . 2001). A crude 
interpolation between the evolutionary tracks yields evolu¬ 
tionary masses of 30.1±1.1 M@ and 22.4±0.8 M 0 for the 
primary and secondary, respectively. This corresponds to a 
mass ratio of 1.34±0.07, in excellent agreement with the 
observed ratio of 1.31±0.02. The orbital inclination corre¬ 
sponding to these evolutionary masses would be 65±2°, and 
this lower value of the inclination would prohibit eclipses. 
A photometric study of HD 93161 is thus necessary in order 
to assess the masses but we might already note that in the 
case of the binary HD 93205 (Paper II), evolutionary masses 
are larger than the true masses. Finally, we caution that the 
luminosity of HD 93161B might be overestimated if the star 
turned out to be another binary system. 


8 CONCLUSIONS 

We have presented a thorough spectroscopic investigation 
of HD 93161, based on medium and high resolution data 
obtained over a period of 7 years. Although the small sep¬ 
aration between the two components of HD 93161 did not 
enable us to get totally uncontaminated spectra, we have 
shown that HD 93161 actually harbours a detached SB2 
system (HD 93161A) and a star displaying RV variations 
(HD 93161B). 

The primary (resp. secondary) star of HD 93161A is most 
probably of spectral type 08V (resp. 09V). Using the RVs 
measured for each component of HD 93161A, we provide the 
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Figure 6. Location of components Al (open circle), A2 (filled 
circle) and B (filled triangle) of HD 93161 in the Ilertzs^rung- 
Russe ll diagram. The evolutionary tracks are from[Schaller_et_alJ 
Il992l) . Isochrones corresponding to ages of 2, 3 and 4 Myr (dotted 
lines) are also drawed. 


first orbital solution for this system. The orbital period is 
8.566±0.004 days and the minimum masses for the primary 
and secondary stars are 22.2±0.6 M 0 and 17.0±0.4 M 0 , re¬ 
spectively. Such large minimum values, compared to typical 
masses for 08V and 09V stars, suggest a high inclination of 
the orbit. It is therefore possible that HD 93161A undergoes 
eclipses. A follow-up study should be undertaken to check 
the photometric behaviour of HD 93161A. 

The visual companion of HD 93161A, HD 93161B, is of type 
06.5V((f)) and is thus slightly hotter than its neighbours. 
HD93161B presents radial velocity variations, but no fully 
obvious signature of a companion could be found in its 
spectrum. The small number of our observations of this 
star did not enable us to investigate the periodicity of the 
velocity changes. Additional observations are required to 
check the binary status of HD 93161B. 

Considering that HD 93161 lies at a distance of 2.5 kpc, we 
found that all three components of HD 93161 are slightly 
evolved off the ZAMS. In the X-ray domain, HD 93161 
appears as a moderately bright X-ray source displaying a 
two-temperature spectrum and no significant variability, 
suggesting an X-ray emission predominantly intrinsic to the 
stellar winds of the individual components. 
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